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1. Introduction 



Light-by- light scattering is a paradigm process in quantum field theory even if it is not 
yet of experimental relevance.^ The calculation of the process at the one-loop level, 
the first non vanishing order, was accomplished quite some time ago [2, 3, 4, 5, 6]. 
Gauge invariance and IR/UV finiteness leads to enormous cancellation mechanisms 
when individual Fcynman graphs arc summed up. As such the study of light-by-light 
scattering is an ideal testing ground for new methods relevant for loop calculations. 

The computation of the higher order corrections to this process has been plagued 
by the difficulty of evaluating two-loop four point functions with four on-shell legs. 
At present it is not possible to include the mass of the particles circulating inside 
the loops. However, the complete set of integrals for planar on-shell two-loop graphs 
with masslcss internal particles is now known [7, 8, 9, 10, 11. 12, 13] together with 
a set of algorithms for reducing the tensor integrals down to the basis set of master 
integrals [14, 15, 16]. This technology has recently been applied to calculate the 
two-loop matrix elements for a wide range of 2 — >■ 2 scattering processes [17, 18, 19, 
20, 21, 22, 23, 1, 24, 25]. 

Recently the two-loop corrections for hght-by-hght scattering were calculated [1] 
by applying the helicity formalism at the two-loop level. The helicity formalism, 
which is strictly defined only in 4-dimensions greatly simplifies the calculation on the 
one-loop level. The generalisation to the two-loop case was done only recently [17, 
23, 1, 24]. In other cases [18, 19, 20, 21, 22, 25] the interference between two-loop and 
tree level amplitudes directly relevant for next-to-ncxt-to-leading order calculations 
has been computed by working with n-dimensional kinematics only. This approach 
is not applicable in the case of four-photon scattering as a tree level operator simply 
does not exist. 

In the present paper we follow another approach to calculate helicity amplitudes 
which in a certain sense interpolates between both methods. Namely by analysing 
the tensorial structure of the amplitude [2], valid to all orders in perturbation the- 
ory, we define projectors to isolate the coefficients of particular tensor structures. 
These projectors are defined in n-dimensions and therefore allow for a completely 
n-dimensional treatment of the computation of the tensor coefficients. The symme- 
try of the process dictates that there are only three independent coefficients. By 
fixing the helicity of the external photons, it is then straightforward to apply helicity 
methods to extract specific helicity amplitudes from the general tensor in terms of 
the three n-dimensional coefficients. The relation between the helicity amplitude and 
the tensor coefficients is independent of the order that the coefficients are computed 
at. It is also independent of whether or not the particles circulating in the loops are 
massive. As an explicit example, we focus on the circularly polarised amplitudes. 
Similar techniques can be applied to obtain linear polarised amplitudes. Note that 

^See the nice review of the current experimental situation in Ref. [1]. 
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in defining the helicity amplitudes the external states are four-dimensional, however, 
the tensor coefficients are fully n-dimensional. 

We apply the method to a fairly general class of Lagrangians containing pho- 
tons, charged and neutral scalars and charged and neutral fermions with prescribed 
couplings. There are several gauge invariant classes and we list these contributions 
separately. We give explicit results for the N — 1 and N — 2 SUSY QED La- 
grangians [26, 27, 28]. Although these unbroken supersymmctric theories have no 
direct phenomenological relevance, they have a theoretical interest of their own and 
have led to new insights in the possibilities of quantum field theory, in particular re- 
garding its divergences. The presence of a supersymmetry leads to important cancel- 
lations between the bosonic and fermionic degrees of freedom. One of the best known 
consequences is the disappearance of quadratic divergences in N — 1 supersymmetric 
theories. These divergences generally appear in renormalization constants and have 
no direct physical consequences. More interesting are cancellations or simplifications 
for physical quantities. These can be analysed on the basis of supersymmetric Ward 
identities [29, 30, 31]. In the case at hand these lead to the prediction that two of the 
independent helicity amplitudes vanish, a fact that we will verify in our calculation. 
The general analysis of supersymmetric Ward identities is a rather involved affair, 
due to the fact that there is no easily implementable supersymmetric regularization. 
The identification of a supersymmetry preserving regulator is beyond the scope of 
this paper. We therefore work in conventional dimensional regularisation which ex- 
plicitly breaks supersymmetry by keeping n — 2 polarizations for internal photons 
while the photino has only 2 degrees of freedom. However, because the quantities 
we calculate are infrared and ultraviolet finite, they are supersymmetric in the four- 
dimensional limit. There may be violations of supersymmetry oi 0{n — 4), but they 
are not relevant in the four-dimensional limit. 

While the simplifications in = 1 supersymmetry are already very interesting, 
stronger results are possible in the case of = 2 supersymmetry. For instance 
it is known that the beta function in A^ = 2 super Yang-Mills theories is fully 
determined at the one-loop level. Since the infinite contributions that determine the 
beta function satisfy such strong non-renormalization theorems, one can also expect 
large cancellations for finite quantities. In order to check for such cancellations 
in the photon-photon scattering amplitude it is necessary to work at the two-loop 
level, since at the one-loop level there is no difference between A^ = 1 and N = 2 
supersymmetry. We will indeed find an interesting pattern of simplifications taking 
place going up in the number of supersymmetries. Similar simplifications for certain 
supersymmetric higher loop amplitudes were already reported in the case of A^ = 4 
Super- Yang-Mills [32] and A^ = 8 Supergravity [33]. 

Our paper is organised as follows. In Sec. 2 we recall the Lagrangian for Super- 
symmetric QED and review the particle content of the theory. The relevant Feynman 
rules are collected in Appendix A. The general tensor structure for light- by-light 
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scattering is detailed in Sec. 3, together with the constraints from transversaUty of 
the external polarisation vectors, Bose symmetry and the gauge Ward identities that 
reduce the number of independent tensor coefficients to three. The general tensor is 
mapped onto circularly polarised helicity amplitudes in Sec. 3.1. We explicitly show 
how each of the three independent helicity amplitudes can be written in terms of the 
three n-dimensional tensor coefficients. In Sec. 3.2 we introduce projection operators 
that can be applied to the full tensor to isolate any of the tensor coefficients. A 
discussion of the Supersymmetric Ward identities is given in Sec. 4.1 while the re- 
sults of an explicit calculation of the one- and two-loop helicity amplitudes are given 
in Sees. 4.2 and 4.3 respectively. The individual contributions to the two-loop am- 
plitudes from the gauge invariant sub-classes of diagrams are listed in Appendix B. 
Finally our findings are summarised in Sec. 5. 



2. The SUSY QED Lagrangian 

N — 1 Supersymmetric QED [26, 27] is an abelian gauge theory containing a vector 
multiplet {A'^, Aq, A") containing the vector photon field and the Majorana photino 
and two chiral multiplets (■?/'£, 0^) and {ip'^, 0^) with charge = — 1 and = +1, 
each consisting of a Wcyl spinor and a scalar field representing the electron and 
scalar electron matter fields. The N = 1 SUSY QED Lagrangian with massless 
matter fields is given by, 

+ \D,<p-L\' + \Dl<i>^\' + mrD,<i/ 

+ V2e {^Pr^kPI - (PWiPr"^ + h.c.) 
-le'{\rLf-\&y, (2.1) 




where 

and where is the gauge covariant derivative and F^,^ is the field strength of the 
photon. This Lagrangian corresponds to the Wess-Zumino gauge where all unphysical 
fields are removed by gauge transformations. 

N — 2 Supersymmetric QED [26] contains in addition a chiral multiplet, a 
Majorana gaugino and a complex scalar photon, (■0^, t/^^", 0°). The N — 2 SUSY 
QED Lagrangian with massless matter fields is given by, 

+V2e (^Pl720Z + (I>rWl^ + h.c.) 
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+V2e (^Pi^0° + 0°*^Pi?^j 
-2e^(|0Zr + |<^+|f |<^°r 

-^e^(l0Zr + l0^.lf + ^e^(l0Zr-|0^.r)'' (2-3) 

where 

Note that the N = 2 supersymmetry leads to a SU{2) symmetry under which the 
photinos and scalar electrons transform nontrivially. By writing them as doublets 

Hi) ' il)' 

the photino-electron-selectron interactions can then be written as 



? + h.c. (2.4) 

72/ 



which is equivalent to the manifestly SU(2) invariant form presented in [28]. 

For completeness, the Feynman rules for these Lagrangians are given in Ap- 
pendix A. 

3. The tensor structure of the four photon amphtude 

We consider the process where all particles are incoming 

7(Pi, Ai) + 7(p2, A2) + 7(P3, A3) + 7(P4, A4) 0, (3.1) 
and photon i carries momentum pi and has helicity Aj. The amplitude has the form 

M = £l,^,£2,;..£3,^3£4,M4A^''^^'''=^^*(Pl,P2,P3,P4), (3.2) 

where the scattering tensor ^^iA'2M3/i4 j^g^ ^j^g following general decomposition 

3 

~ \ 3\3i y -f^ji t'di ~ nji y ■'^ji ■'^j2 ~ .71.72 y ■f^.71 ■f^.72 

+ ■ o^2''3 V^-^ V^^ + P^ • q'"'^" V^^ V^'' + P*' • Q^''^" V^^ V'"'\ 

^ 3132 y i^3l ^32 ^ 3132 y ^Ji ^^J2 ' 3132 ^ ^3i ^32 J 
3 

+ X] ^3l32j33'lPji PhPh Pjl- (3-3) 
jlj2j3j4=l 
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Terms such as eA'iM2At3Ai4 g^j-^ forbidden on the grounds of parity invariance. The 
functions Ai, and Cijki are functions of the Mandelstam variables, s = (pi+P2)^, 
t = (p2 + PsY, u = (pi + and the spacetime dimension n. They also depend 
on the mass of the electron and any other particles that may be involved in the 
scattering process. This decomposition is valid for arbitrary loop order since it is 
based solely on the Lorentz structure of the external particles. Altogether, there are 
138 coefficients. However, many coefficients are irrelevant since they drop out when 
contracted with the photon polarisation vectors due to the transversality condition, 

ej-p, = 0, forjG {1,2,3,4}. (3.4) 

This reduces the number of coefficients to 57. 

Bose symmetry of the external photons means that these coefficients arc not in- 
dependent. Requiring invariance under the exchange of any index pair (j, /Uj) where 
the first labels the external vector and the second is the Lorentz index of the corre- 
sponding polarisation vector reduces the number of independent functions to 11. 

The number of independent functions is further reduced by gauge symmetry. 
The related Ward identities read in an obvious notation, 

^pi£2e3£4 _ ^£1P2£3£4 _ ^£1£2P3£4 _ ^£1£2£3P4 _ Q (3 5) 

Applying the gauge symmetry reduces the number of independent functions to three 
which we take to be, 

Ai{s,t,u), Bl-^{s,t,u), C2iu{s,t,u). (3.6) 

Once these functions are known, the full tensor can be reconstructed. 

3.1 Helicity amplitudes 

It is often convenient to express the amplitude in terms of the helicity structure of 
the scattering process. This can straightforwardly be achieved using the Lorentz 
structure of the tensor. In principle there are 16 helicity amplitudes - two polarisa- 
tions for each photon. However, parity, time-reversal and Bose symmetry reduce this 
number to four which is further reduced to three by crossing symmetry. Note that 
the number of independent helicity amplitudes matches the number of independent 
functions describing the tensor. We take 

M+++-, M++-, (3.7) 

to be a basis from which the other helicity amplitudes can be reconstructed. 

By specifying polarisation vectors for the external photons and applying them to 
the full tensor, we can derive the helicity amplitudes in terms of the three independent 
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functions Ai, Bl^ and C2111. We find that, up to overall phases, 

M++++ = Ai{s, t, u) + Ai{t, u, s) + Ai{u, t, s) 

--B\^{s,t,u) - ^-^-±B\^{s,u,t) - -Bl^{t,s,u) 

^ '-B\{t,u,s) + ^-—^B\{u,s,t) - -^-—^Bl,{u,t,s) 



— 2t)M^ , , , stit — u)^ , , t^U^ /, X 
+ ^2111^, M + ^^2111 s, u, t) + — C2111 s, u) 

st^ ttt^ Sli^ 

~'-^C2iii{t,u,s) + — C2iii(w,s,t) + —C2iii{u,t,s), {3.1 
ou OS or 

= -— 5i\(s, -u) - u, t) + -^Bliit, s, u) 

s(t — u) , , , „i , X -sf— t + m) ^1 , 

s) + ^i?u(«, s, + ^-^^Bliu, t, s) 

s(t — u)u^ , , , Stit — u)^ , t^M^ /, X 

H — 62111 (s,t,M) 62111 (s,M,t) - — 6211l(t, s,m) 

81 8k 8s 

St^ ^ , , tu^ ^ , , SM^ 



— 62111 (i,ii,s) - — 62111 (ii,s,i) + — 62111 (ii,t,s), (3.9) 
M 8s 81 



and 



— = Ai(s, t, u) + Ai{t, u, s) + Ai{u, t, s) 

, X t{2s-\-t). , X „1 / 

--B\^{s,t,u) - ^-^-±B\^{s,u,t) - -BUt,s,u) 

sit — u) . , , (2s — t)t . , , s(s — 2t) . , , 
^ '-B\^{t, u, s) + ^—-^Bl,{u, s, t) - -^-—^Bl,{u, t, s) 



s(s — 2t)u^ , , , stit — u)^ , t^u ^ , , 

+ 62111 s, t, u) + ^62111 s, u, t) + — 62111 t, s, M 

8r 8it 8s 

si?' tip' sii} 

--^62111(^,14,5) + — 62111 (M,S,t) + — 62111 (M,t,s). (3.10) 
QU 8s 8t 

Amplitudes for linearly polarised light can also be straightforwardly obtained from 
the general tensor and are given in terms of Ai, B\^ and 62111. 

Note that in deriving the helicity amplitudes, we have made no assumptions 
about how the functions Ai, B\^ and 62111 are calculated. In the conventional 
approach to computing loop helicity amplitudes one has to define a description to deal 
with scalar products between loop momenta and polarisation vectors. Furthermore, 
we have made no assumptions about the masses of particles in the loops. 

3.2 Projection operators 

To calculate the independent functions Ai, B\y and 62111 it is convenient to define 
projection operators that (a) isolate the function and (b) saturate the free Lorentz 
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indices. This latter point is important for practical calculations since it allows for the 
cancellation of reducible scalar products between loop momenta and external vectors 
in the Feynman diagram integrals. These reducible scalar products can be expressed 
in terms of inverse propagators and many cancellations already happen at this level. 
For example at the one-loop level, no tensor box integral needs to be evaluated since 
any loop momentum appearing in the numerator will immediately cancel one of the 
propagators. This reduces the complexity of the one-loop calculation enormously. 
For higher-loop or multi-leg applications such simplifications are highly appreciated. 

To invert Eq. (3.3) and isolate the independent functions, it is useful to define 
the following tensors. 



j=l k=l 

k=l 



(3.11) 
(3.12) 



where H — G ^ is the inverse of the 3 by 3 Gram matrix defined by the momenta 
Pi-i P2, Ps- In terms of Mandelstam variables, these matrices are given by 



G 



f s u\ 
s t 
u t 



-4 



( 1/s + l/u 
1/s 
1/u 



1/s 



(3.13) 



1/u 

1/s + l/t 1/t 
1/t 1/t + 1/u I 

P is a projector onto the {n — 3) dimensional subspace perpendicular to the 3 di- 
mensional space spanned by the vectors pi, Ps, Pi- 'T^j,v is the dual vector to 
relative to this 3 dimensional space. One may easily check the following relations for 
the objects V^v and 7?.^, 



0, 

-P^, = tr{V) 
V^Wj = 0, 



tr{V • P) = n - 3, 



(3.14) 
(3.15) 
(3.16) 
(3.17) 
(3.18) 
(3.19) 



To determine the tensor coefficients, Aj, we define the following linear operators 



ii(M) = 
i3(A^) = 



1 


(n — l){n 


-3) 


1 




(n-l)(n 


-3) 


1 





(n-l)(n-3)^''^''^^''^'^^-^^ 



(3.20) 
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One finds by direct calculation that 

tr{V) 



(MM)\ 

A2{M) 

[mm)) 



(n- l)(n-3) 



(tr{V) 1 1 \ 

1 tr{V) 1 
I, 1 1 tr{V)) 



A2{s,t,u) 
\A3{s,t,u) J 



(3.21) 



Staying in n-dimensions and inverting the system of equations yields, 



^ Ai{s,t,u)\ 

A2is,t,u) 
\As{s,t,u) J 



(n-4) 



[{n-2) -1 -1 \ (Ai{M)\ 

A^iM) 

[mm) J 



\ 



-1 
-1 



(n-2) -1 
-1 (n-2)y 



(3.22) 



We see that the right hand side appears to introduce an additional factor of l/(n — 4). 
On the other hand, the sum 

Ai{s, t, u) + ^2(5, t, u) + ^3(5, t, u) = Ai{M) + Ai{M) + Ai{M) 

is free of spurious poles. We will discuss how these apparent extra poles cancel each 
other directly in the helicity amplitudes below. 

It is straightforward to define projectors for the coefficients B^^, 



BliiM) 



1 



(n-3) 
2 

" (n - 4) 



Hki {{n - 2)ii(M) - A^iM) - A,iM)) , (3.23) 



where the last term is proportional to Ai. The other i?"'s are found by permuting the 
Lorentz indices in (3.23) and subtracting the respective A^s. Acting with Eq. (3.23) 
on the general tensor one finds the desired property, 

Bti{M)^B^,{s,t,u), 

by using relations (3.14-3.19). 

In the same spirit one can construct projectors for the coefficients Cjkim, 

CjMmiM) = nj^,nk^,ni^,nm;.,M'''^'^'^' 
-2 



Hj.BUM) + HjiBl^iM) + Hj^BUM) 

%{M) + H,^B% 



+HuBUM) + HkmBUM) + Hi^B%{M) 



-4 



^'"rS ((^ - ^)MM) - A^iM) - As{M)) 



(n-4) 

.^^^ ((n - 2)A2{M) - As{M) - A,{M 

■^^^ ((n - 2)AsiM) - A,{M) - A^iM 
[n — 4:) ^ 



(3.24) 
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where the last three terms are proportional to Ai, A2 and respectively and which 
again satisfies, 

Cijki{M) = C\jkiis,t,u). 

The linear operators acting on the amplitude defined in Eqs. (3.20), (3.23) and 
(3.24) are sufficient to determine any of the coefficients on the right hand side of 
the general tensor decomposition of the amplitude in Eq. (3.3). The linear nature 
ensures that one can apply these projectors on a graph by graph basis. 

All of the tensor coefficients are hnearly related to the basis set of Eq. (3.6) by 
gauge invariance and permutation symmetry. An important and non-trivial check of 
a full diagrammatic calculation is to evaluate each of the tensor coefficients and to 
verify the relations amongst them imposed by the Ward identities. 

It remains to be shown that the poles present in Eq. (3.22) cancel. This is seen 
by focussing on the dangerous Aj terms only. For C2111 one finds that the l/(n — 4) 
pole drops out directly, 

C211l{s,t,u) = C2Ul{M) 

^ ... + 4HuHi2 \Ai{s, t, u) + A2{s, t, u) + A^is, t, u) 



For Bl^ one has. 



Bl,{s,t,u) = Bl,{M) = ... + —A,{s,t,u). (3.25) 



Using the Bose symmetry relations, 

^i(s, t, u) = Ai{s, u, t) , A2{s, t, u) = Ai{u, t, s) , ^3(5, t, u) = Ai{t, s, u) , 

a short calculation shows that in Eqs. (3.8), (3.9) and (3.10) only the sum of the A/s 
is present so that the helicity amplitudes are free of spurious poles. This ensures that 
the projector method can be applied on a graph by graph basis without the need to 
expand one order higher in e. 



4. Light-by-light scattering in SUSY QED 
4.1 The SUSY Ward Identity 

Supersymmetric amplitudes are expected to satisfy the Supersymmetric Ward iden- 
tity [29, 30, 31], 

= (0| [g,4...at] |0), (4.1) 

where Q is the Supersymmetry generator that satisfies Q|0) = and a\ arc the 
creation operators for particles in the initial state. Application of Eq. (4.1) to states 
with three photons and a photino yields the following identities, 

= 0, (4.2) 
= 0. (4.3) 
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At tree level these identities are trivially satisfied due to the abelian nature of the 
photon. Beyond tree-level, couplings with the matter multiplet give rise to non- 
trivial interactions, that we investigate in the subsequent subsections. As mentioned 
in the introduction, we work in conventional dimensional regularisation and treat all 
momenta and Lorentz indices as n-dimensional. This breaks supersymmetry because 
the fermionic and bosonic degrees of freedom for the gauge multiplet are no longer 
equivalent. We therefore expect that there will be non-trivial corrections to the 
SUSY Ward identity that vanish in the 4-dimensional limit. 

For convenience we decompose the helicity amplitudes perturbatively as 

Aa.A3A. = O? {m^SmX^M + f -^ilAaA, + ^(«')) ■ (4-4) 

4.2 One-loop SUSY QED helicity amplitudes 

Because of the abelian nature of QED, there is no quartic photon coupling at tree 
level. At one-loop a four-point contribution is generated by the interaction of the 
photon with the matter multiplet. The one-loop amplitude can be decomposed 
according to gauge invariant subsets of graphs. There are two such groups at one- 
loop so that 

=Ma).5 + ^(i),i^^ (4.5) 

where the dependence on the helicities has been suppressed. These two contributions, 
M^^^'^ and Al^^^'^, denote the scalar electron and electron loops respectively. In the 
Standard Model, only the electron loops, At*^^)'^, contribute. We find 

M^lY++= 8 + 0(6), 
M2i^+ = -8 + 0(6), 

M'iJ-= 8 + 0(e), 

= 4 - 2 ((A: - Y)2 + TT^) + 4 ((A: - Yf + TT^ - 2{X - Y)) ^ 

S 

+ |m -l-O(e), 

M^IY— = -4 - 4 ((a: - Yf + tt" - 2{X - Y) 

+ |ii^t|+0(e), (4.6) 
where 

X = log (—), Y = log (—) . (4.7) 



s / \ s 

Combining all graphs together, the one-loop light-by-light scattering amplitudes 
in SUSY QED are rather compact and are given by, 

^ 
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Figure 1: Typical diagrams involving the photino that are related by charge conjugation. 
The dashed line represents the scalar electron, while the overlaid solid and wavy lines 
represent the photino. Similar graphs can be drawn where the dashed line represents the 
scalar photon and the overlaid solid and wavy lines represent the electron. 



M 



(1) 



= 



= 4s{n - 4) (Box%s, t) + Box%s, u)) - 4s(n - 2)Box''(t, u) 

= -4(^{X-Yf + ^) + 0{e), (4.8) 



where Box^(s,t) is the (infrared and ultraviolet finite) one-loop box graph in 6 — 2e 
dimensions. At one-loop there is no contribution from the gauge multiplet and this 
is therefore the result for both = 1 and = 2 SUSY QED. 

The fact that A^Ji++ and both vanish is directly attributable to the 

SUSY Ward identity. The zeroes for Ai++++ and occur at the level of 

the master integrals, i.e. to all orders in n — 4. This is perhaps not surprising 
since although dimensional regularisation breaks the supcrsymmctry for the gauge 
multiplet, the photon and photino arc not present as internal particles in any of the 
one-loop graphs. At this order, we are not sensitive to the SUSY breaking engendered 
by dimensional regularisation. 

At one-loop we expect that amplitudes contain terms of weight 0, 1 and 2 (count- 
ing logarithms and tt as weight 1, squares of logarithms as weight 2 and constants 
(excluding tt) as weight 0). By inspection of Eq. (4.6), we see that all such terms 
are present in the fermion and scalar contributions. However, the supersymmetric 
cancellation is such that only terms of weight 2 remain in Eq. (4.8), all weight and 
weight 1 contributions are eliminated. We also note that terms proportional to ratios 
of the kinematic scales t^/s^ in the individual contributions (4.6) completely cancel 
in the supersymmetric result of Eq. (4.8). 

4.3 Two-loop SUSY QED helicity amplitudes 

Unlike the one-loop case, the chiral photino-electron-selectron coupling is present and 
it is necessary to address the issue of how to treat 75. Because of the parity invariance 
of the process, the scattering amplitude contains no antisymmetric part, indicating 
that the 75 contributions drop out. We can see this by considering a typical pair of 
diagrams with a photino exchange that are related by charge conjugation as shown 
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in Fig. 1. Up to overall factors, the contribution from the first diagram is given by 

f dPk Tr(PR/PLAglAlg2Al2g3Al23g4A) 



where ki_j — k -\- Pi-\- . . . -\- pj and 



(4.10) 



Similarly the second graph is 



' y (27r)" (27r)" {k-')\k^Y{k^A)\k2MYe-\k -if ^ ' ^ 

RelabeUing k — > —A;, I ^ —I and using charge conjugation to reverse the trace, we 
find that, 

d'^k d'^i Tr(Pi/PR/c£:iyfcie2Ai2£^3Ai23^4A) 



(27r)« (27r)- (A;2)2(A;i)2(A;i2)2(A;i23)2£2(A; - £)2 
so that, up to common factors 

/i + /2 - WPl + Pl/^Pr] ■ ■ ■ = ^ [/ - 75/75] • • • • (4.13) 

As expected, traces with single 75 factors drop out entirely. Since the amplitudes are 
finite, we therefore use an anti-commuting 75 such that, 

^ [/ - 75/75] = /. (4.14) 
Similar arguments apply to the chiral couplings of the scalar photon with the electron. 
4.3.1 N = 1 SUSY QED 

As in the one-loop case, it is convenient to break the amplitude up according to the 
different gauge-invariant pieces so that 

= + >^(2)-^ + ;t<(2),P + j^my^ (4.15) 

where the dependence on the helicities has been suppressed. At two-loops there are 
contributions from photino exchange, A^*-^-*'^, and the four scalar vertex, A^'^^),^'^ 
well as graphs where the electron or scalar electron couple directly to the photons, 
and Al^^),^. Altogether there are 1902 Feynman graphs which we generate 
using QGRAF[34]. Figures 2, 3, 4, 5 and 6 show the types of Feynman graphs 
relevant for the various contributions. Tadpole graphs and self-energies of external 
legs are not shown since they vanish in dimensional regularisation. We note that all 
of the possible diagrams are planar. 
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A 



B 



C 



Figure 2: Graphs relevant for the electron, photino, scalar electron and scalar photon 
contributions M.^"^^'^, M.^"^^'^, M.^'^^'^ and TW^^)'^. The solid lines represent particles from 
the matter multiplets, i.e. the electron or scalar electron while the internal wavy lines 
denote particles from the gauge multiplets, the photon, the photino or the scalar photon. 




D E F 




G H 

Figure 3: Graphs relevant for the photino and scalar electron contributions M.^^^'^ and 
J\^4i'^)'S , The solid lines represent the scalar electron while the wavy lines denote either the 
photon or photino. 





J 



K 





O 



Q 



Figure 4: Graphs relevant for the scalar electron contribution A^(^)''^. The solid lines 
represent the scalar electron while the wavy lines denote the photon. 
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R S T 

Figure 5: Graphs relevant for the scalar four-point interaction contribution M^'^^'^ . The 
solid lines represent the scalar electron while the wavy lines denotes the photon. 




U V W 

Figure 6: Onc-particlc reducible graphs relevant for the electron, scalar electron and scalar 
photon contributions Ai^'^^'^ , M^'^^'^ and Ai^"^^'^ . The solid lines represent particles from 
the matter multiplet, i.e. the electron or scalar electron while the wavy lines denote the 
photon or scalar photon. Graphs V and W are only part of the scalar electron contribution. 

Explicit evaluation of the two-loop graphs yields rather lengthy results for the 
individual gauge invariant subsets, and we list them in Appendix B. Combining 
the various components according to Eq. (4.15) gives the following two- loop = 1 
SUSY QED helicity amplitudes, 

= + 0(6), 
MiV = + O(e), 

MS- = +(| - 16 U,{y) + 8^Y + 2^ Us{y) + 8 ^^3(0;) 

+8 (^XY-Y + X^U2iy)-l6C3 + ^ - ^ (gY + 11 + 2XY^ n'' 

+ - X"^ + 6 X^Y"^ + AXY^ + 12 XY - AX^ + SY^yj 

+ 1-16 LU{z) -S[Y-x] Lisiy) -S[Y-x] Us{x) - 8 XYLi^iy) 



+ ( -^^^-SXYjU2{x) + ^7T''-^ [Y + 2-X]{Y-X]n' 
1 



+ ^ l^AX^ + UX^Y + 12XY^ -UX'^Y -AXY^ + Y^-3X^ 
J J u 

+4"^(5 (v--A-).^ + i(v-Ar)^ 
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+ ( 16 Li3(x) + 8 X U2{y) + 8X U2{x) 




u^t\+0{e), (4.16) 



where we have used the standard polylogarithm identities [35] to express our results 
in terms of a basis set of constants (where is the Riemann Zeta function, ^2 = '^''^/S, 
Cs = 1.202056 . . .), logarithms and polylogarithms Li„(w) defined by 



Un{w) ^ / — Li„_i(t) forn = 2,3,4 
Jo t 

rw rjf 

U2{w) = - -log(l-t). (4.17) 
Jo t 



with arguments x, 1 — x and {x — l)/x, where 

x^ — , y = — = l-x, z^--^ . 4.18 

S S t X 

In the physical region s > and t,u <0, our basis set of functions are all real. 

As expected from the SUSY Ward identities, Eqs. (4.2) and (4.3), two of the 
helicity amplitudes vanish in the n — > 4 limit. However, the identities are violated by 
terms of (9(e) which can be traced back to the SUSY breaking nature of dimensional 
regularisation. 

The SUSY Ward identity does not require that vanishes, but, as in 

the one-loop case, there arc still significant cancellations and the full amplitude 
is somewhat more compact than the individual contributions (which are listed in 
Appendix B): In fact, at two-loops we expect that amplitudes contain terms up to 
weight 4 (counting Li„, In" and Cn ^ weight n, ln"*Li„ as weight n + m and so 
on). Each of the individual contributions for listed in Appendix B, and 

(2) F 

in particular the fermion contribution A4\.+ that corresponds to the Standard 

Model {N — SUSY), demonstrates that all possible weights, 0, ... ,4, arc present. 
However, we see that the = 1 SUSY amplitude of Eq. (4.16) contains terms 
of weight 2, 3 and 4 and that, as in the one-loop case, terms of weight and 1 
cancel. We also note that the individual contributions contain terms proportional to 
the dimensionless ratios, f^/v^s^, i'^/s^i i/u (and t m), while in the combination 
selected by A^ = 1 SUSY, (4.16), the terms proportional to t^/u^s^ drop out. 

4.3.2 A^ = 2 SUSY QED 

For N — 2 SUSY, there are additional contributions from the gaugino, the scalar 
photon and the modified quartic scalar electron interactions. In terms of the different 



15 



gauge-invariant pieces we find that 



where again the dependence on the hehcities has been suppressed. Here, A^(^)'"^ 
denotes the 144 graphs involving the scalar photon while the factors of 2 and 3 mul- 
tiplying the photino and four-point scalar electron contributions respectively reflect 
the aditional photino and the modified scalar electron interactions of the N — 2 
SUSY theory. 

Combining the individual gauge invariant contributions listed in Appendix B 
according to Eq. (4.19), we find the N ^ 2 SUSY QED helicity amplitudes are 
rather compact and are given by, 

AlS++ = + O(e), 
M?|+_ = O + 0(e), 

MS— =1-16 Lu{y) + 8YUs{x)+8Y Us{y) + —n^--XY7r^ 





+i7r| 1^16 Lk{x) - -Ytt' - -Y' \ Y -3X 

+ i^u^t^ + 0{e). (4.20) 

As expected, two of the helicity amphtudes vanish as n — > 4 due to the SUSY Ward 
identity. The remaning non-trivial helicity amplitude is considerably simpler than 
that obtained in either pure QED [1] (see the fermion loop contributions in Ap- 
pendix B) or the N — 1 SUSY QED hehcity amplitudes of Eq. (4.16). In particular, 
we note that only terms of weight 4 remain, the contributions of weight 2 and 3 (that 
were present in the = 1 case) have cancelled. Furthermore, all terms depending 
on the ratios of kinematic scales have dropped out. 
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5. Summary 

We have demonstrated that the method based on n-dimensional projections is able 
to generate hehcity amphtudes in an efficient way. We have been able to confirm the 
previous results for photon-photon scattering via a charged fermion loop obtained 
by helicity methods in the high energy limit where the fermion mass can be ignored. 
The method can, in principle, be used for more complicated processes, involving 
massive particles in the loop, non-abelian fields and/or more external vertices. The 
method is constructed such that one can use standard n-dimensional Lorentz covari- 
ant reduction techniques for tensor-integrals not only for squared matrix elements, 
but also for helicity amplitudes. 

As an application we studied photon-photon scattering in the theoretically in- 
teresting cases of = 1 and N — 2 supersymmetric QED. Because the process is 
both ultraviolet and infrared convergent, the results are explicitly supersymmetric 
in n = 4, which is sufficient to draw conclusions regarding SUSY cancellations. We 
have not addressed the question of whether the method can be used to define super- 
symmetric amplitudes also at higher order in (n — 4). However, we expect that the 
situation here will be similar to that found using other techniques, where explicit 
terms proportional to {n — 4) have to be added in order to preserve the SUSY Ward 
identities at all orders in (n — 4). 

The results we found have an interesting pattern of simplification as one increases 
the number of supersymmetries. In normal QED, corresponding to A" = super- 
symmetries, one finds three independent amplitudes. For supersymmetric theories, 
two of these hehcity amphtudes vanish. This is well understood on the basis of the 
supersymmetric Ward identity. An intricate pattern of cancellations occurs in the 

remaining (non-trivial) amplitude, , as A^ is increased, both in terms of the 

dimensionless ratios of kinematic scales and in the weights of the functions present 
in the amplitude. 

First, considering the one-loop level, one finds in normal (A^ = 0) QED terms 
proportional to the dimensionless ratios t^/s^ and j together with weights of 
logarithm of 0,1, 2 (see Eq. (4.6). In the supersymmetric case (4.8), there are no 
dimensionless ratios and there is a uniform weight of 2 for the logarithms. 





one-loop 


two-loop 


A^ = 


0,1,2 


0,1,2,3,4 


N ^\ 


2 


2,3,4 


N ^2 


2 


4 



Table 1: Weights of terms contribution to M.+-\ at one- and two- loops. The N = 

result refers to the fermion contribution of Eq. (B.3) while values given for N = 1 and 
N = 2 are extracted from Eqs. (4.16) and (4.20). 
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t U t 11 


/4 ,,4 


N^O 


V 








V 






N = 2 









Table 2: Ratios of kinematic scales appearing in the two-loop amplitude M.+-\ for 

different amounts of supersymmetry. 

At the two- loop level, the pattern is similar, see Tables 1 and 2. Normal (A^ = 0) 
QED contains powers of t/s up to 4 and all weights of logarithm from to 4. The 
N = 1 case shows some simplifications, powers of t/s exist only up to squares and 
the logarithmic terms have weights 2, 3, 4. Finally in the N = 2 theory there are no 
powers of t/s in the amplitudes and the logarithms are uniformly of weight 4. We 
see therefore, that increasing the number of supersymmetries reduces the complexity 
of the amplitudes and that the N — 2 theory is maximally simplified. A conjecture 
would be that this pattern persists in higher orders of perturbation theory, where 
one would expect no powers and a weight of two times the level of loops in the 
logarithms. This pattern might be explainable by performing similar calculations in 
superspace. 
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A. The SUSY QED Feynman Rules 

The Feynman rules for the SUSY QED Lagrangians of Eqs. (2.1) and (2.3) with all 
momenta incoming are given by, 




'Pi 



P2 





2ie^ i = j 

-ie^ i^j,N = l 



-ie^/2PR 
-ieV2PL 
+iey/2PL 
-ieV2PR 



e 

e 
e 
e 



71 

72 



7 

7 

(l>R: 7 = 7i 

7 = 72 



(A.l) 




f -ieV2PL 
-ieV2PR 
+ieV2PR 
-ieV2PL 



7 
7 



7i 



e 

e = 0S, 7 = 7i 
e = 0S, 7 = 72 




-ieV^PR $° = 



where the photon, electron, scalar electron, photino and scalar photon are denoted 
by wavy, solid, dashed, overlaid solid and wavy and overlaid dashed and wavy lines 
respectively. 
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B. Two- loop contributions to the helicity amplitudes 

In this appendix we list the individual contributions to the two-loop helicity am- 
plitudes. There are five separate gauge invarariant contributions, loops with scalar 
particles coupling to photons, electron loops, photino exchange graphs, diagrams 
with the four scalar vertex and graphs involving the scalar photon. We find that 
individual two-loop contributions are given by, 

MtiU = -6 

= -12 

= +24 
= -12 

A^S+_ = -4X2-^+ (2Y + XY -SX^^"^ - 4x(2 + 3X^ ^ 

+ (^-2 71^-2 (2Y -2X -2XY -SX"^ + 3Y'^yj^ 
+i7r|-8 1^1 + 3x1 - + 4 ( -1 + 2X^ - 8 ( - 2 X + Y^ ^ - 8 X ^ 

-2X^-^- (-X^ + 2y) -2x(2 + 3Xy- 
+ (^- TT^ - (2Y -2X -2XY -SX"^ + 3Y^y^ 



+i7r<^-4 1 + 3X - + 2 X-1 -4 -2X + Y\—-AX 



\ J U \ J \ J u^s^ 

+ (3 7r2 + 3 \ 2Y -2X -2XY -^X"^ + ?>Y'^\\^— 



s2 



+i7r|l2 fl + 3xl--2 f5X-3) +12 { - 2X + y\K + 12 X ^ 



U \ I \ I v?s'^ 



+ < -u ^ t 
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iX'4^-8x(x + l]- 



+ 1^16 hu{x) -s(y -2^ Li3(y) - 8 Li3(x) Y -s(x-Y + XY^ U2{x) 



16C3-§7r^ + ^ (^-3r-3-X2 + 2Xrj7r2 
1 



+^ I - 24X2 -3- 18X^^2 _48j^^gy3^24XF 



» 

+ ^ - 16 hu{z) - 32 Li4(|/) - 32 U^^x) 

+8 |^-2 + y + 3xjLi3(|/) + 8 |^F + 3X + 2j hh{x) + S U2{y) XY 



+ f- ^7r2 + 8(^-2y-2X + Xy^^ Li2(x) + ^ vr^ 



+- ( -6xy + 4y + y2 + 3-2X + x2 1 



3 
4 
3 



Qx^ + x^ + ex^y-gx^y^ -4^3 y + 6xy2 -6y- 2x^ + 2 



+6x + 6xy))^ 

+i7r| -8 1^1 + 2 ^ - 8 X 



+ - 16 Lialy) - 8 X Li2(?/) - 8 Li2(x) X + - tt^ y 

+^ |^6y2-i2-6xy-9xy2 + y3jj 

+ 1^32 Lisd/) + 32 Li3(a;) + 32 y Lis^/) + 32 [y - ij Li2(a;) 
^-2 + 5y + 3X^ TT^ 

|^-9x2y + 6y2-27xy2 + i2xy + 6y-i2x-6x2 + y^ + 3x^ 

-2^2 ^ 



+ i^Q U^{y) + s(y-X^ U2{x) - 16 Ca - ^5 + 6 yj tt^ 



23 



+^ l^ - 6Y^ + 12XY + AY^ - 3 - 12Yjj -4^^ 

+ 1^-16 Lu{z) - 16 Lu{y) - 16 Lhix) +8 (^-1 + 2X^ Usiy) 

+8 (^1 + 2X^ Lk{x) + - ^ TT^ - 8 + j Li^ix) + ^ TT^ 

(y"^ + 2XY -Y -X -S + X^^TT^ (^5X^ -6XY^ -12Y -WX^ + Y^ 
-18X^ -AXY^ -20X^Y + 6X^Y^ + 12XY + 12X + 10Y^ + 30X^Y] 



2 



+i7ri -A -+ (t TT^ + 8 ( - XY + Y -1 + yA\ -A X ^ 



+ (^16 Li3(y) + 16 Lisix) + 16 YU2{y) + 16 - 1^ U2(x) 

(^-l + 2Y + 2X^7r^ 
+^ (l2 XY^ - SY + 6 X^Y - 3Y^ - 2 X^ + 6 X + 3 X^ - 6 Xy\\ ^ 



+ \u^t\ (B.3) 



A^S^_ = 6X2 4^ + 4x(4X + 3)- 
U^S'' \ J u 

+ (^-32 Lu{x) + 16 Us{y) Y + 16 (y - ?j Us{x) + 8 (y - X^ U2{y) 

47 2 / \ 
+ 16 U2{x)YX + 16 C3 + ^ TT^- 3 ieXY -2X'^ -6X -bW'^ 



— ( -i5x2y^-i8y^ + i2xy-6-36y-4x^y + y^ + 4y' 



+ 1^16 LU{z) + 48 Lu{y) + 48 Lu(x) - 8 y- 3 + 2Y + AX jLk{y) 

-8 ^2 y + 4 X + 3^ Us(x) - 16 U2(y) X Y 

+ (^^ 7r^-8 (^-3Y + 2XY-3Xyj U2{x) - y tt^ 

(llY + 2Y^ -16XY -5X + 9 + 2X^^77^ 
+^ ( -18Y + 3X^ + 15XY^ -27X^ -7X^ + Y^ -30X'^Y^ -AXY^ 



24 



-12X^Y + 1SXY + 1SX + 7Y^ + 21X^yJ^ ^ 
I 3 + 8X1 - + 12 X ^ 



+ 1^32 Li3(a;) + 16 X Li2(y) + 16 U2{x) X F + 1 j tt^ 

+ 1^-48 Lisd/) - 48 Li3(a;) - 48 F U^^y) - 48 - 1 j Lis (a;) 
+^ ^-3 + 8y + 4X^ 7r^ + ^ |^2r^ + 9y + 9y^ - 18X 

-12^2 y -42xy2 + 4x=^- + i8xyjj ^| 

= -6 X^ 4^ - 12 X f 2 X + lU 

+ (^2 hu{y) - 16 X Li3(x) - 16 X U^{y) - g vr^ - ^ - 1 + 2 X^ j tt^ 
+ ^ 1^-3 -4X^^-18^2- 3X2^2 ^r^-36rjj 

+ 1^-48 Li4(y) - 48 Li4(a;) + 24 |^1 + F + xj U^{x) + 24 + X - 1 j U^{y) 
-24 I^F + xj Li2(a;) + ^ 7r^ + ^ ||l3F + 3X2 + 9 - 7X + 3^^ - 12XFj tt^ 

^- i2X=^r-36y + 36X- 10X^^ + 3x^-54x2 + 30x2 y + 3y^ 

+10y^ - 12Xy=^ + 42Xy2 + 36Xy _ I8x2y2j j ^ 



+i7r|-12 ^1 + 4xj ^ - 12 X 



t 



4 



+ ( -32 Li3(y) + - X7r2 + - ( - 9 + - 3X2 y - 9 F 



+ (48 hh{x) + 48 Li3(y) - 48 U2{x) - 4 ( y + X - ij 7r2 



i 

-4 |^6X y + x^- 3x2- 3xy2-6X- 3x2 y + 3y2 + 3y + y'^j J — 



t2 
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MjiL = - 1 (B.4) 

The contribution from fermion exchange A1^^^'^ has previously been computed in 
Ref. [1] and we find complete agreement with the results presented there once the 
different definitions of the helicity amplitudes are taken into account. 



26 



